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Advanced  Simulation  of  Metal  Clusters 

Summary  of  effort 

In  an  effort  to  develop  efficient  0(N)  first-principles  calculation  algorithms,  we 
devised  the  approach  using  wavelets  as  basis  functions  for  electronic  structure 
calculation.  The  theory  of  wavelets  allows  one  to  apply  a  multi-scale  (multiresolution) 
analysis  to  problems  that  exhibit  widely  varying  length  scales.  Furthermore,  the  dual 
localization  property  of  the  wavelet  bases  is  useful  for  improving  the  existing  0(N ) 
methods  that  are  yet  based  solely  on  the  real-space  locality.  Our  application  of  this  idea  to 
the  density-functional  molecular  dynamics  method  leads  to  a  robust  algorithm,  which 
holds  potential  to  extend  the  applicability  of  the  current  first-principles  methods  to 
systems  an  order  of  magnitude  larger.  On  the  other  hand,  we  have  developed  an  approach 
to  the  development  of  reliable  “pair-functional”  interatomic  potentials  for  monoatomic 
metals  based  on  fitting  to  a  large  set  of  experimental  and  ab  initio  data.  In  connection 
with  the  construction  of  improved  and  accurate  semi-empirical  models,  we  have 
investigated  the  simulation  tools  such  as  hyper-Molecular  Dynamics,  generalized 
simulated  annealing  for  global  optimization,  finite-difference  method  for  calculating  non¬ 
linear  static  third-order  polarizabilities. 

Accomplishments/New  Findings 

(i)  Parallel  and  Distributed  Molecular-Dynamics  Simulations 
Parallel  Algorithm  and  Simulations 

We  have  implemented  a  parallel  and  distributive  algorithm  for  classical  molecular 
dynamics  using  many-body  “pair-functional”  potentials,  and  tested  a  parallel  version  of 
ab  initio  molecular  dynamics  algorithm.  The  former  parallel  molecular  dynamics 
algorithm  was  tested  for  systems  consisting  of  1 06- 1 07  atoms,  and  was  employed  for  the 
simulation  of  the  epitaxial  island  growth  on  hexagonally  reconstructed  Au(100).  The 
study  for  Au/Au(100)  reveals  that  the  stable  islands  of  rectangular  shape  are  hexagonally 
reconstructed  in  conformity  with  the  patterns  of  the  reconstructed  Au(100)  surface,  and 
suggests  the  “magic”  width  of  the  islands  in  agreement  with  experimental  observations. 
The  associated  study  on  transition  states  and  adatom  diffusions  indicate  that  the 
experimentally  observed  strong  anisotropic  effect  is  attributed  to  the  long-range  exchange 
diffusion. 

Glassy  Surfaces  of  Metallic  Nanocrystals 

The  vibrational  and  elastic  properties  of  metallic  nanocrystals  are  studied  using 
semi-empirical  potentials  of  gold  and  lead.  Our  results  on  the  density  of  vibrational  states 
show  that  the  experimentally  observed  linear-frequency  dependence  is  attributed  to  the 
glass-like  surface  atomes,  as  a  result  of  a  thorough  geometry  optimization.  The  extracted 
elastic  constants  indicate  a  decrease  of  Young’s  moduli  and  an  increase  of  the  shear 
moduli  with  the  increase  of  the  nanocrystal  size.  The  changes  involved  in  the  elastic 
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constants  can  be  attributed  to  well-defined  surgface  layers,  independent  of  the 
nanostructure  size.  The  experimentally  observed  enhancement  of  low-  and  high- 
frequency  modes  are  found  to  be  spatially  localized  in  two  shells,  a  soft  surface  shell  and 
a  hard  transition  shell. 

Nonlinear  Optical  Properties 

We  have  devised  a  method  for  calculating  valence  electron  contributions  to  the 
static  molecular  third-order  polarizabilities.  The  method  utilized  is  based  on  the  finite- 
field  approach  coupled  with  semi-empirical  polarization  calculations  on  all-valence 
electrons.  The  static  third-order  polarizabilities  of  carbon-cage  fullerenes  are  analyzed  in 
terms  of  three  properties,  from  a  geometric  point  of  view:  symmetry;  aromaticity;  and 
size.  The  application  to  large  fullerenes  shows  that  the  static  linear  polarizability  depends 
almost  exclusively  on  surface  area  while  the  third-order  polarizability  is  affected  by  a 
combination  of  number  of  aromatic  rings,  length  and  group  order,  in  decreasing 
importance. 

Hyper  Molecular  Dynamics 

Recently,  based  on  the  transition  state  theory,  Voter  has  provided  a  method  to 
accelerate  the  molecular  dynamics  (hyper-MD),  which  opens  a  window  to  simulate 
atomic  dynamics  for  a  microsecond  or  longer.  In  this  scheme,  a  bias  potential  raises  the 
potential  energy  except  for  the  saddle  points  of  the  potential  energy  surface.  The 
dynamics  on  the  biased  potential  surface  leads  to  accelerated  evolution  from  one  potential 
minimum  to  another.  We  have  explored  the  possibility  of  using  a  local  bias  potential  to 
speed  up  the  hyper-MD.  This  approach  is  shown  to  be  a  good  approximation  to  the  real 
system.  So  far,  the  hyper-MD  works  only  for  interatomic  potentials  with  continuous  third 
order  derivatives. 

(ii)  Multiresolution  Analysis  for  Electronic  Structures 

Orthonormal  Wavelet  Bases  for  Electronic  Structure  Calculations 

We  have  devised  an  approach  utilizing  compactly  supported,  orthonormal  wavelet 
bases  for  ab  initio  molecular  dynamics  (Car-Parrinello)  algorithm.  A  wavelet  selection 
scheme  is  developed  and  tested  for  various  atomic  and  molecular  systems.  The  method 
shows  systematic  convergence  with  the  increase  of  wavelet-selected  grid  size,  along  with 
improvement  on  compression  rates.  This  method  yields  an  optimal  adaptive  grid  for  self- 
consistent  electronic  structure  calculations,  and  offers  a  realistic  approach  for  the  study  of 
metal  clusters. 

The  application  involves  the  study  of  atomic  and  electronic  shell  structures  for  the 
experimentally  characterized  AI77.  The  onion-like  AI77  structure  can  be  described  as  a 
stable  All 3  inner  core  covered  by  a  two-layer  atomic  shell.  The  stability  of  AI77  is 
confirmed  by  the  structural  optimization  and  electronic  structure  calculations.  The 
interaction  between  AI77  and  its  ligand  is  found  to  be  ionic-like.  Our  results  on  the 
electronic  structure  provide  important  information  on  the  electronic  shell  structures  in 
large  metallic  clusters. 
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Highly  Accurate  Electronic  Structure  Calculations 

The  structure  and  energies  of  the  (1x2)  and  (1x3)  missing-row  reconstructed 
transition  metal  (110)  surfaces  are  studied  using  first-principles  density-functional  theory 
with  local-density  and  generalized  gradient  approximations.  A  detailed  analysis  of  the 
multilayer  relaxation  shows  that  the  heavy  topmost-layer  distortions  are  closely  coupled 
with  the  change  of  kinetic  energy  of  electrons.  The  first-principles  results  shed  lights  on 
the  nature  of  roughening  and  deconstruction  transitions  of  Au(l  10). 
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